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Abstract The symmetric simple exclusion process is one of the simplest out- 
of-equilibrium systems for which the steady state is known. Its large deviation 
functional of the density has been computed in the past both by microscopic and 
macroscopic approaches. Here we obtain the leading finite size correction to this 
large deviation functional. The result is compared to the similar corrections for 
equilibrium systems. 
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1 Introduction 

Over recent years there has been a growing interest in understanding the fluctua- 
tions and the large deviations of the density of systems in a non equilibrium steady 
state |[Tl2l3T4l5l^^ . In such steady 

states, the generic situation is that the correlation range of density fluctuations 
extends through the whole system 0221123 Il24ll231l26ll271l28l|29l and the large de- 
viation functional of the density is non local l2l[3lffll5ll^l8ll9l [T0l[TT1[l^[T3ll . This 
contrasts with systems (with short range interactions and far from a critical point) 
at equilibrium, where the range of correlations is microscopic and the large devi- 
ation functional is local. 
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Two main approaches have been followed recently to study these large de- 
viations: for some microscopic models such as exclusion processes the steady 
state measure is known exactly ll8ll^ [TO[TTl[T2l[T3l and finding the large deviation 
functional is a matter of computing large scale properties (very much like when 
one tries to calculate the free energy in equilibrium systems starting from the 
Gibbs measure). Obviously this microscopic approach is limited to cases where 
the steady state is exactly known. The other approach is the macroscopic fluctu- 
ation theory J2][3]|U[5][6]EI) f° r diffusive systems which calculates the large devi- 
ation functional by identifying the optimal path followed by the system to gener- 
ate a given deviation. In systems at equilibrium, time reversal symmetry gives a 
simple relation between this path and the relaxation path starting from the same 
deviation, and so identifying this optimal path is easy. On the contrary in non equi- 
librium steady states this time reversal symmetry does not hold and the approach 
is limited to cases where the equations giving this optimal path can be solved. 

One motivation to study the large deviation functional of the density is that it 
generalizes the notion of free energy to non equilibrium states 171 1131 . As analytic 
expressions of these large deviation functionals are usually hard to obtain, they 
are known so far for a rather limited number of models. The one dimensional 
symmetric simple exclusion process (SSEP) was one of the first models [368 
1911131 for which an explicit expression could be derived which showed the non 
local character of this large deviation functional. The goal of the present work is 
to obtain the leading finite size corrections to this large deviation functional and 
to compare it with the corrections one typically finds in equilibrium systems. 

The SSEP describes a lattice of L sites in which each site i is either occupied 
by a single particle or empty B3II34| [35 36 37 38] . Each particle independently 
attempts to jump to its right neighboring site, and to its left neighboring site with 
rate one. It succeeds if the target site is empty; otherwise nothing happens. At the 
boundary sites, 1 and L, particles are added or removed: a particle is added to 
site 1, when the site is empty, at rate a, and removed, when the site is occupied, 
at rate y; similarly particles are added to site L at rate 8 and removed at rate /3. 
These injection and removal rates at the boundaries correspond to the left and right 
boundaries being in contact with reservoirs at densities 

a 8 

(to justify fl} it is easy to check using detailed balance that, when one forbids the 
exchanges of particles at site L by setting f} = 8 = 0, the steady state measure is a 
Bernoulli measure where all the sites are occupied with probability p a . Similarly 
one can check that when the contacts between site 1 and the left reservoir are 
broken, the system equilibrates at density pi,.) The main advantage of the SSEP is 
that its steady state measure is known [23 30 8,9 1 for arbitrary a,/3,y, 5 and L. 

Here we try to determine the generating function of the density, which is sim- 
ply the Legendre transform of the large deviation functional. Let P{n\, ■ ■ ■ ,ni) be 
the steady state measure of a one dimensional lattice gas on a lattice of L sites, 
where «, > is the number of particles on site i (for the SSEP the only possible 
values are «, = or 1 but in the more general case discussed in section 4 the occu- 
pation numbers n, will be arbitrary). We want to calculate the following generating 
function Zt(h\ , ■ ■ ■ ,hi) (which, in the equilibrium case, is nothing but a partition 
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function in a varying field) 



Z L (/ii,--- ,h£) = £ ex P kf"!'"'"!) (2) 

{«,,-»,} V=i / 

where hi depends on the site i. Let us define Gi(h\,-- ■ hi) as 

G L {h u ---h L )=\ogZ L {h u --- ,h L ) ■ (3) 

We would like to obtain an expression of Gi(h\ , • • • , hi) for a slowly varying field, 
that is when h\ is of the form 

h >= H (i-k) (4) 

where X is large (the reason for the shift of — 1 /2X in © is simply to make sites 
1 and L play symmetric roles). This choice for the dependence of hi allowes one 
to test density deviations which vary on a large length scale X, which might be 
different from the system size L 

In the following we will consider the case where the two lengths X and L are 
large (compared with the lattice spacing) but comparable 

L = Xy (5) 

with y of order 1 . For the symmetric exclusion process in contact at site i = 1 and 
at site i = L with two reservoirs at densities p a and pi, it is known [8 , 9 EH ( see 
for example eq. (80,81) of lfT3l ) that, in the steady state, 

G L (h u ---h L )~X min frfi(log(l+F(i)( e "W-l))-logfj-^) 

{F(x)}Jo V V Pb-PaJ 

(6) 

where the minimum is over all the monotone functions F(x) which satisfy F(0) — 
p a and F(l) = pi,. From © it is easy to see that the equation satisfied by the 
optimal F is 

„ F'(x) 2 (e H ^ - 1) 

F " {X) = l-F(x)+F(x)e H W (7) 

Note that the non local character of the functional © comes from the fact that the 
optimal F(x) depends on all values of H(z) for the whole range < z < y = j. 

On the other hand for a system at equilibrium (with short range interactions) 
one expects lp7l[T3l that 

G L (h lr --h L )~X [ y dxg(H(x)) (8) 
Jo 

where g(h) = lim/ j _ i . 00 Gi(h, - ■ ■ , h) /L is the extensive part of the free energy in a 
constant h. This is obviously a local functional of H(x). 
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The main result presented in the present work is that the leading correction to 
© is 

G L (h h -h L ) ~ A j^x(log(l +F(x)(e ff W - 1)) -log(gg)) 

-«log (£g) -blog - ilog( 9 (0)) (9) 



^Pb-Pa J B \Pb-Pa) 2 

where T 7 is the solution of (0, the parameters a and b are defined as in 191 1271 

a = — !— ; b=—^— (10) 
a+Y P+8 

and (p(x) is the solution of the linear differential equation 

♦'w — (£$)'*<*> <"> 

which satisfies the boundary conditions (p(y) = and (p'(y) = — 1/y. 

This can be compared to the case of a systems at equilibrium where the form 
of the leading orders of Gl is 

G L {h , • • -h L ) ~ A f g{H{x)) dx +A Mt (H(0)) + A^ ht (H (y)) 
Jo 

[ H'(0) B le{t (H(0))+H' (y) 5 ri s ht (H (y)) 



X 



./o 



C{H{x)) H'{xf dx 



where A left (/!),A ri s ht (/z),fi left (/j),B ri s ht (/j) and C(A) are defined in (146152153b . We 
see that the leading correction (i.e. the term of order in A) is also non local (0 in 
the out of equilibrium SSEP whereas it corresponds to boundary contributions A left 
and A nght in the equilibrium case (112b . At the next order (the order 1 /A), which 
we did not study in the non-equilibrium case, one can notice in ( fl2l ) an integral 
containing the gradient term H'(x) 2 characteristic of the Ginzburg-Landau theory. 

The paper is organized as follows. In section 2, we do a direct perturbative cal- 
culation when the hj's are small and we check that the expansion agrees with the 
prediction (0 for the SSEP. In section 3, we present the derivation of (0 for ar- 
bitrary hj's. In section 4, we discuss how (112b can be derived in the equilibrium 
case. 



2 Perturbations for small h t 

In this section we present the straighforward calculation of Gl from the knowledge 
of the correlation functions of the density in the steady state. From the definition 
©, one can relate the expansion of Gl in powers, of the hi's to the steady state 
correlations. For example to second order in the hj's one has 



G L (h u —h L ) =£fti(«i)+£y«n?)-(«i) 2 ) + £^ (ninj) c + 0(h 3 ) (13) 

i i '<j 
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For the SSEP, it is known that in the steady state H23II271 



M = (nf) = PaiL + b -J ) lf; a ~ l \ (14) 
L + a + b — 1 

, > (p fl -Pfc) 2 (i+a-i)(L+fc-y) 

For large L and A (keeping their ratio constant as in ©, when the hi have 
the form © the various sums in (1131 can be computed by the Euler Mac Laurin 
formulae: 



e f/(^-|)=/;V(^-^ i^£izm +0(£ 



Lf'(Le)-/'(0)] , _ 4 , 



£2 1< § < /( /e "l) g (-' £ "f) = f^ dx s{y)dy - y o L£ f(x)g(x)dx+o( £ 2 ) 
and one gets 

JXm) = A [ y p(x)H(x)dx+(p a - p b ) [" y - 2ay - 2 * + 2 2aX + 2bx H(x)dx (15) 
; Jo Jo 2v z 

Lr2; \ i P-/ \„2/ \ j . I \ /" v y-2flV-2x+2flX+2to 2 
Af<«i> =A/ p(x)H z (x)dx+(p a -p b ) -= H 2 xdx 
Jo Jo 2y z 

2>?(„,.) 2 = X /V(,)fl*(,)«fo+ ( Pa - Pb ) f y - 2ay - 2X + 2aX + 2bX p{x)H 2 {x)dx 

j Jo Jo y 4 

where p (x) is the steady state profile 

p (x) = p°<y 

and 

«,« j \niriji c = a - pb) j 
■ (Pa-Pfc) 

( 1 - 2a) (.y - x) (.V - z) + (4a + 4fc - 6)jr(y - z) + (1 - 26)xz 

We have checked that these expressions coincide with (|9]l at second order in H (x). 
For example at first order in H(x) the solutions of Q and dTTT ) are 



Y/1,/2 =-X(Pa-p b f [ y dx P ' dz X -^^ L H(x)H{z) 

77; Jo Jx y 

y P x(y-x)H{x) 2 dx+ P dx P dzH(x)H(z)x 

J0 Jo Jx 



\ p a {y-x) + p b x 2 

F(X) = + {pa - Pb) 



P X -^H(z)dz- P Z -^H(z)dz 
Jx y 3 Jo y 3 

<p(x) = (y-x)/y+ [Pa ~ Ph) P(2z-y-x)H(z)dz 
y Jx 



{P a-Pb) 

y 2 Jx 

and inserting these expressions into (|9) one gets (115b 
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3 Derivation of the main result © 

Our approach to obtain (|9| consists in choosing h, piecewise constant: takes n 
possible values H\,- ■ H n in n consecutive boxes. As in each of these boxes, h t 
is constant we will use the expression (116117b for a single box with a constant h 
which is much easier to obtain. Then we will use an additivity formula ( 171b to go 
from the expression for one box to the expression for n boxes. Finally we will take 
the limit n — > <*> to establish ©. 



3.1 A single box 

Using the matrix ansatz (see the Appendix lAl. one gets, for large L, the following 
expression for Zi(h,h, •••/?) by dividing ( 169b by i 



Z L (h,--,h) 



where 



{pa-pb) L+a+b ^ L - a ~ b 

\+ Pa {e h -\)Y(\+p b {e h -\)f 
l+p B («*-l) 



^ = ^r-r lo g 



1 °l+p b {e h -l) 



(16) 



(17) 



Remark: Let us check that these expressions agree with the claim (|9) in the intro- 
duction: one has by solving Q and dl lb for a constant h 



F(x) = 



(1+ Pa (e h - 1))W> (l+p b (e h - !))*/> _ l 



and 



<pw = 



y-x 



from which it follows that 



log[l + (e A -l)F(*)]-log 



yF'{x) 

Pb - Pa 



= log 



Pa-Pb 

Mo 



y(e h -l) 



\+ Pa {e h -\) 



F'(y)= l+ . p f- 1) i°s fl+Pb{eh - 1) 



y{e*-\) 



l+p a (e"-l) 



and by replacing into (0 one finds an expression equivalent to (116b obtained by 
the direct calculation. This shows that (O is valid in the case of a constant hi. 
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3.2 Several boxes: the prediction © 



Let us now come to the case of several large boxes with a constant h, in each box. 
We will first write down the expressions predicted by the claim (O. Then we will 
see in the next subsection that these expressions coincide with those obtained by 
a direct microscopic calculation. 

For piecewise constant H(x), with 



H(x) = H m for x m -\ < x < x„ 



with 



Xq — , X m — x m ~ \ ~\-ym 
the solution of Q in the interval x m _i < x < x m is 

1 



x n =y 



(18) 
(19) 



Fx = 



(1 + (e H - - (1 + (e H - - \)F m 



i - 1 



(20) 

where F m = F(x m ). Writing that F'(x) is continuous (i.e. F'{x m ) = F'{x m ) + ) at 
all the x m 's leads to the n — 1 equations that these F m 's should satisfy 



y m (e H --l) 



log 



\ + {e H »< -\)F m 
I + (<?"• -l)F m - 



l + ( e ",„+i _i)F m 
y m +l{e H »'+i -1) 



log 



l + ( e H»,+i -l)F m+1 



(21) 



Equations ( 120121 1 ) fully determine the solution of (0 for a piecewise constant 
H(x). 

To solve the equation dl 11 1 for <p(x), one can first notice that the discontinuity 
of the <p'(x) at x = x m is 

9 (x ffl ) + - <p (x,„) = — — f(x m ) . 

f [X m ) 

Everywhere else the function (p(x) is piecewise linear. These jumps of <p'{x) and 
the fact that (p(y) = and <p'(y) = —1/y determine the function (p(x) everywhere: 
in the interval x m < x < x m +\ one gets 

<p(x) = -\y-x + £ (y-Xmi)(x mi -x)W„ 1 (22) 
y m\>m 

+ £ {y -Xnti ) {Xm } - *m 2 ) (*m 2 ~ x)W mi W mi 
my >m 2 >m 

+ £ Cy-* mi )(x mi -Xm 2 ){Xm 2 -Xm 3 )(Xm 3 -x)W mi W m2 W m3 + ...] 



mi >m 2 >ni-i>m 

with W m defined by 



F 1 [Xm) 



Using ( 1201 ) and (1211) . one can show that 



W m 



■log 



IF, 



m— 1 



■log 



1)F 



m+1 



' 1 + (e Hm 
l + (gtf,„ + i_l)F m 



(23) 



(24) 
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In summary in the case of several large boxes the claim ((9) leads to 

Z L (hi,-- A) = e Gi(,,1 '- /!i) ~^exp[A<*f(p fl ,p fo )] (25) 

with 

l + (e H >"-l)F m 



^ = "^ miog L(P,-p")(^-D iog (i 



m-l 



(26) 



and 

« _ / Pb_Pa_Y ( Pb^PaV 1_ 

V^W \yF'(y)J <p(0)V2 U/J 
with F(x) and <p(x) given by (120122b and the F m 's solutions of d2TT) . 



3.3 Several boxes: the microscopic approach 

Let us now see how the microscopic calculation for the single box can be gener- 
alized to the case of several boxes and leads to expressions equivalent to ( 126127b . 
We consider the case of several large boxes with a constant h, in each box 

hi = H m for Li + ---+L m -i<i<Li + ---+L m (28) 

uses the additivity formula d71b and the saddle point method. 
Let us define 

^_ (p a ,a\(e Hl D + E) Ll \p,b) 
Zl(Pfl ' Pj - (p a ,a\(D + E^\p,b) 

and for i > 2 

{p,\-b\{e H -D + E) L '\p' 1 b) 



Zi{p,p') = 



(p,l-b\(D + E) L i\p',b) 

which are the generating functions for each box (see Appendix |A). (z\ is special 
simply because in Z2 , • ■ ■ Z n the parameter a has been replaced by 1 — b). Then using 
the additivity formula (171b derived in Appendix lAl one gets for Zi(h\,--- ,h£) 
when 

L\ +L2 H L n = L 

and the hi are of the form (1281 

Z ^-^) = r(L + a + b) X (29) 

[dpi f dp n -x (p a - Pb) L+a+b Zl(pa,pl) Z2(Pl,P2) • • • Z„(pn-l,Pb) 
] lilt" '1 2iK {pa-pl) Ll+a+b (pl-p2) L2 + 1 ---(pn-l-pb) L " + 1 

where the integral contours verify pt < \p tt -i \ < • • • < |Pi | < Pa- So far ( |29l l is 
exact for arbitrary H m 's and L m 's. The virtue of (171b is that it relates the properties 
of the whole system of those of the n subsystems. 
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When the lengths L,„ of the boxes become large, if we define the y m 's by 

L m = A y m 

one knows ( fT6b from the single box calculation that 

Zm (p, P ') ~2Upy)^ m[Am(p ' p ' )+log(p ~ p,)1 » 

with 

A m (p,p') = -log 



1 , l+pTe^-l) 
■log 



and for i > 2 



Bi(p,p') = 



B«(P,P') 



1 & l+p'(e"<»-l) 

(p -p') n + fo e (" +z, ) A i(P-P') 
(l+p(e"i -l)) a (l+p'{e H i-l)) b 

(p-p')<A( p ' p,) 



(1 +p(e«' - l)) 1 " fo (l +p'{e»i - l)) h 
Then using the saddle point method in (129b one finds that 



^* e A[.c/*+ylog(p a -p 6 )] 



(30) 
(31) 

(32) 

(33) 



where 



.s/* = min \y 1A1 (p fl , n ) +y 2 A 2 (n , r 2 ) + ... + y„A n {r n -i,p b )]-ylogy+Yyi log y t , 
{'•») S 

(34) 



(Pfl-Pfo) fl+6 Ji^ 2 5i(P fl ,n) yIS 2 (ri,r 2 ) ylB n (r n ^,p b ) ^ +m ^-i/ 2 



and 4 is a tridiagonal matrix 



n - r 2 



r n -i -Pb 



(det[4])" 



4 = 



/(/, Vi \ 
V! c/ 2 v 2 o 
... 
0.. 
0. 

V o o o o y„_ 2 £/„_! / 







(35) 



(36) 



with 



U„ 



(P^m+l)] 



dp 2 

a 2 [y m+ lA m+ i(p,p')] 



dpdp' 



(37) 



(38) 



p=r m ,p'=r m+ i 



(Note that in 04b one takes the minimum and not the maximum because the inte- 
gration contours are perpendicular to the real axis, and a maximum over r m along 
a contour becomes a minimum when r m varies along the real axis). 
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The saddle point values r\, ■ ■ ■ r n -\ (i.e. those which achieve the minimum in 
d34b ) satisfy the n — 1 equations 

d [y m A m (r m _i,p) +y m+ iA m+l (p,r m+l )] _ Q 

d P p=r a ~ 

These saddle point equations turn out to be the same equations as those satisfied 
( 121b satisfied by the F m 's. Therefore one has 

r m =F m (39) 

This already allows one to verify, using (130134b . that the term proportional to A in 
d25l l and ( 1331 is the same. 



One can also show by a direct computation that the U m 's and the V m 's defined 
in ( 137138b can be expressed in terms of the function F(x) given in d20b . 

( \ 1 \ 1 F"(jc m )(-) -F"(x m )W 

\y m y m +\JF'{x m y F'{x m y 

Vm = ~pr~, \ j— i / / r ■ (41) 

y m +\F '{x m )F '{x m+ i} 

It is then easy to see that one can rewrite U m as 

1 1 \ 1 W m 

U m 



Jm Jm+1 / F f (x m y F 

with W m given in ( 1231 ). Then the determinant of the matrix A defined in 
be computed 

det [ 4 ] = ^77—! iFTT 7f7 x ( 42 ) 

yi---y„[F'{xi)---F'(x n -i)\ z 

y ~\~ ( y x m ) x m W m -h ( y x m j ) {Xffi j x m2 )x m2 W m j W m2 

m m\>mi 
m\ >m2>m-$ 

and using the fact (see (13012011 ) that 

exp[A (F { F m )}= l + = 1 + (^"-1)^-1 (43) 

ymF {Xm) ymF (x m _j) 

one gets for 38* 

b 



Pb-Pa\ Pb-Pa 



yF'(0) J \yF'{y 

I -j- y Xm ^ ~ Xm ^ \y _|_ ^ (j ~ x mi ) { x m\ ~ x m 2 ^/ 
m ^ mi >?«2 y 



-1/2 



This expression coincides with the expected expression (127b . Therefore this sub- 
section has established the validity of (O in the case of several boxes. 
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3.4 A large number of boxes 



Let us now try to take the large n limit of the above result. We consider that we 
have n boxes, of equal length L/n = Xy/n, and that the field H m in the m-th box is 
given by 

'my y 
n 2n/ 

where H (x) is a smoothly varying function. For simplicity we choose the boxes of 
equal lengths. Therefore one has 



H m — H 



x m 



my 



One then need to solve the equations d2TT l satisfied by the F m 's. For large n one 
can show by a direct calculation that the solution of these equations is given by 



<ny\ 
n J 



O 



1 



(44) 



where F(x) is the solution of 0. (Note that from now on, F(x) is the solution 
of when H{x) is a smoothly varying function. This solution F{x) is not identi- 
cal to ( 120b which was obtained for a piecewise H(x). The difference is at the origin 
of the correction of order 0{nT 2 ) in d44b . This difference will lead to negligible 
terms anyway.) One has from (1301 



m I "m— 1 s r m ) 



log -+ log 
n 



l+F(x m )(e H M-l) \ 
-F'(x m ) J 

y F{x m )H\x m )e H ^ 
In 1 + (e»M _ \)F{x m ) 

which can be rewritten using the fact that F(x) is solution of (0 



O ^ 



A m (F m - h F m ) = -log-+log 
n 



l+F{x m )(e 



H(x„ 



-F'{x m ) ) 
\+F{x){e H ^ ~\) 



O - o 



Using then the Euler McLaurin formula to perform the sum (1341 . one finds that 



gf(pa,p b ) = -nlog-+ / logf 
n Jo \ 



l+F(x)(e 



Hx) 



1) 



-F'(x) 



dx + O 



(45) 



and this leads ( 133134b to the term proportionnal to A in (0. 



One can also obtain the large n estimate of W m 

w - y ( F " i -"Y 
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Then by defining W(x) by 

F"(x)\ 



W(x) 

one can see that 

s _ {P a -p b ) a+h x 
(-F'(0))«(-F'(1))* 



F'(x) 



1 /■! r-.,_7Ur_^v \-V2 



1+ r^l W (x)dx + I dx [ d z {y -^ Z - X)X W{z)W{x) 
Jo y Jo Jx y 



Then if one <p(x) is solution of 

<p"(x)=W(x)<p(x) 
with (p(y) —0 and (p'(y) = — 1/y one has 



yW = i-^+r& (y " z)(z " 3r) w( z ) 
y Jx y 



and one finds 



as claimed in ((9). 



-jw) \^W ] (9(0)) 



4 Equilibrium case 

Let us consider a one dimensional lattice gas on a lattice of L sites, where each 
site i is occupied by an integer «, > number of particles. We assume that the 
interactions are short range and that at equilibrium the system is homogeneous in 
the bulk with correlation functions decaying exponentially fast with the distance. 
We would like to obtain an expression of G]_{h.\,--- ,hi) defined in ( 1213b for a 
slowly varying field of the form (|4) when both X and L are much larger than the 
range <jj of the correlations between the occupation numbers n,. 



4. 1 The constant field case 



Let us discuss first the case of a constant field (hi = h). If g(h) is the extensive part 
of the free energy Gl 

G L (A, •••/») 
8(h) = lim j 
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one expects that in the large L limit 

G L (h) = Lg(h) +A left (h) +A^ ht (h) + o(exp[-L/^ (h)]\ (46) 

where A left (/j) and A nght (h) represent the contributions of the left and right bound- 
aries respectively and ^ (h) is the correlation length (in presence of the constant 
field h). The form ( 146b can be easily understood by the transfer matrix method, in 
particular exp(— 1/^ (h)) is the ratio of the two largest eigenvalues of the transfer 
matrix. These two contributions A left and A nght are not necessarily equal as they 
may differ if one imposes different boundary conditions at the two ends. 
In a constant field, one can also define the average density (n,) at site i by 

d\ogZ L 

(n ' ) = ^r (/v " /j) 

and the pair correlation function 

d 2 \ogZ L 

In the large L limit, far from the boundaries (i.e. for i 3> % (h) and L — i^>£, (/?)), 
the average density («, ) has a limit independent of 

(ni)^g'(h), (47) 

the pair correlation function («;H_,-) C depends only on the distance j — i 

(n,-n y ) c ->■ cj-iQi) (48) 

and one has 

g"(h) = £ c k (h) . (49) 

jt=-oo 

On the other hand, close to the left or to the right boundary, i.e. as long as 
i ~ E, (h) or L — i ~ (h) these quantities keep in general a dependence on even 
in the large L limit. For example 

( ni )-g'{h)^aY\h) ; (n L -i)-g'(h)^af sU (h) (50) 
4.2 The slowly varying field case 

Now for a slowly varying field of the form (|4), when X » ^ (h) (more precisely 
X » max, ^ (hi)) one expects i71 [T3ll that to leading order 

G L (h h ---h L )~X [ } g(H(x))dx 
Jo 

where y is defined in (0. This can be easily understood by cutting the system of 
length L into many subsystems of size Xdx much larger than ^ but much smaller 
than X. In each of these subsystems the field hi is essentially constant, and the free 
energies of these subsystems can simply be added. 
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As explained in Appendix |Bj the leading corrections to this formula when L 
and A are much larger than the correlation length t, (h) is 

G L (h,---h L ) ~ A fg(H(x)) dx + A Mt (H(0)) + A n ^( H (y)) 
■Jo 



H'(0) B Mt (H(0)) +H' (y) B^ ht {H (y)) 



Jo 



C(H(x)) H'{x) 2 dx 



0l F 



which is the result announced in d 1 2b where 



i=\ 



g'(h) 
24 



I 

/=o 



and 



(51) 

af ght (h) 
(52) 

(53) 



k>\ 



Remark: for the same system on ring with periodic boundary conditions, imply- 
ing in particular that H (y+x) =H(x), the boundary terms disappear and one gets 

G L {h u ---h L )~X Pg(H(x))dx+ r P ' C{H{x))H'{xf dx + --- (54) 
Jo A Jo 



Note that at order 4- the integral of H'{x) 2 in d5TT l and in (l54l is nothing but the 
square of the gradient of the Ginsburg Landau theory. 



5 Conclusion 

In this paper we have obtained the first correction (|9) to the large deviation func- 
tional of the density for the non equilibrium steady state of the SSEP and compared 
it with the corresponding term for equilibrium systems dT2l . Like in the equilib- 
rium case, this first correction does not depend on the system size. On the other 
hand in the non-equilibrium case (0 the correction has a non-local character, very 
much like the leading term. Our derivation is based on the knowledge of the steady 
state measure, as given by the matrix ansatz. 

An interesting question would be to try to recover our result by the macro- 
scopic approach: in the macroscopic fluctuation theory [2 3 6], the large deviation 
functional of the density is given by the contribution of the optimal trajectory in 
the space of all the time dependent density profiles which produces a given devi- 
ation starting from the steady state profile. A natural question would be to try to 
calculate the correction by integrating over all the profiles in the neighborhood of 
this optimal profile. Such an approach was successful in understanding the first 
corrections to the large deviation function of the current H 3 9 11401 . and it would 
be of course interesting to see whether it works as well for the deviations of the 
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density. If this is the case, one could try to determine similar corrections for other 
models such as generalizations of the SSEP H12lll5ll4"n . 

Recently, it has been noticed that the large deviation functional could exhibit 
phase transitions [42 s 43 [44]. Whether the corrections calculated here would be- 
come singular at such phase transitions is another question one could try to inves- 
tigate. 

It has also been shown that the SSEP, in a non equilibriuml steady state, could 
be mapped by a non local change of variables onto a system at equilibrium 1145 1 . It 
would be interesting to know whether this transformation could be sued to confirm 
our prediction (|9) and establish a connection with (fT2l . 

Acknowledgements We would like to thank Vincent Hakim for helpful discussions. 



A Additivity formulae 

It is known 13013 11131 and has been used in several previous works 18191 that the steady state 
measure of the SSEP with injection an d remo val rates a , /? , y, 8, as defined in the introduction, 
can be calculated by the matrix ansatz. [30 31]. The probability of any microscopic configuration 
{ni ni} (with iij = or 1) is given by 



(W\X l X 2 ...X L \V) 
P{{nu-n L })= mD + E)L{v} (55) 



where each matrix X, depends on the occupation «, of site i 

X i =n i D + (l-n i )E (56) 
and the matrices D.E and the vectors \V), {W\ satisfy the following algebraic rules 

DE-ED=D+E 

{W\{aE-yD) = (W\ (57) 
(pD-8E)\V) = \V) . 

Given these algebraic rules, one can define a family of left and right eigenvectors (p a ,a\ and 
\Pb,b) by 

(p a ,a\(p a E-(l-p a )D) =a{p a ,a\ (58) 

((1 - p b )D - Pb E) \p b , b)=b \p b ,b) . (59) 
The vectors (W| and |V) which appear in J55l57t are examples of such eigenvectors 

(W\ = (p a ,a\ , \V) = \p b ,b) (60) 



when p a = a/(a + y), p b = 8/{8 + j5) and a = l/(a + y) , b = 1/(5 + j8) as in dlllOt . 
Then for < b < 1 and p a > p b , one can prove the following key additivity formula 
{p*,a\XiX z \p b ,b) rdp (p a -p b ) a+b (p a ,a\Xi\p,b) (p, 1 -b\X 2 \p b ,b) 



{p a ,a\p b ,b) J 2ni (p a - p)"+*(p - p b ) (p a ,a\p,b) (p,l-b\p b ,b) 



(61) 



where X[ and X 2 are arbitrary polynomials of D's and E's and the contour is such that p b <\p\< 

Pa- 
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Proof of <6I> : Let us first derive of the following identity j9] 

(p a ,a\(D + E) L \p b ,b) r(a + b + L) 



{p a ,a\p b ,b) r(a + b) (p a -p b ) L 



(62) 



To do so one can notice that in the steady state of the SSEP, as defined in the introduction, the 
average occupations satisfy 

a-(a + y){n\) = (m -n 2 ) = ■ ■ • (n; - n i+1 ) = ■■■ = (p + 8)(n L )-8 

These L equations which express simply that in the steady state the current is conserved, can be 
solved. From the solution J14t one can see that 

(n,'-n,- + i) = ■ 



L+a+b- 1 

On the other hand using the matrix representation (155157160) one has 

_ ■ _ (p a ,a\(D + E) i - l (DE-ED)(D + E) L - i - , \p b ,b} 
[Hi (p a ,a\(D+E)Lp b ,b) 

= (p a ,a\(D + E) L - l \p b ,b) 
{p a ,a\(D + E) L p b ,b) 

These two identities give the recursion 

(p a ,a\(D + E) L - l \p b ,b) = (p a -p b ) 
(p a ,a\{D + E) L \p b ,b) L + a + b-1 

which establishes the veracity of l l62l . 

Now to prove 1 161) (as in fT3l ) one can first notice that the discussion can be limited to X\ 
and Xj of the form 

Xi = [p a E-(l-p a )D]P' [D + E}'» 

X 2 = [D + E}Pi [(l-p b )D-p b E]<K 

as any polynomial in £>'s and £"s can be written as a sum of such terms (this is because D and 
E are linear functions of the operators A and B defined by A = D + E and B = p a E — (1 — p a )D 
and that AB — BA = A. Thus word made up of A's and B's can be ordered as a sum of terms of 
the form B'^A'^ or A^-B^-) Then the left hand side of (61} becomes 

{p U) a\XlX 2 \p b ,b) = ^ m (p a ,a\{D + Er + i"-\p b ,b) 
{p a ,a\p b Jj) (p a ,a\p b ,b) 

while the right hand side of 1 161) becomes 

API b <n I H_ (P«-Pi>y +h (PaM(D + EY"\P,b) ( P A-b\(D + E)P2\p b ,b) 
J 2ni(p a -pY+ b (p-p b ) (p a ,a\P,b) (p,l-b\p b ,b) 

and the equality of l l63t and d64l follows from the expression 462b and the Cauchy theorem. This 
completes the derivation of <611 . 

First consequence of (ST\ : 

It is possible to show directly from the algebra ( 155157) that 

(p a .a\D\p b ,b) _ bp a +ap b _ (p a ,a\E\p b ,b) _ b(l - p a ) +a{\ - p b ) 
(p a ,a\p b ,b) p a -pb ' {Pa-a\pb,b) Pa~ Pb 
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which becomes by replacing p a by p and a by 1 — b 



(pA-b\D\p b ,b) 



= b- 



Pb 



(p,\-b\p b ,b) p-p b 
Therefore d61) becomes after integration 



{p,\-b\E\p b ,b) = _ bj \-p b 



(p,l-b\p b ,b) 



{p a ,a\X D\p b ,b) =b (pg,a\X \p b ,b} | d_ 
(p a ,a\p b ,b) (p a ,a\p b ,b) dp 



{p a ,a\X E\p b ,b} = b {p a M\X \p b ,b) 

(p fl ,a|Pfe>M {p a ,a\p b ,b} h dp 



Pa-p b Y +b 
Pa- P ) 

Pa ~ Pb 



P-Pb 



(p a ,a\Xo\p,b) 



P=Pi, 



Pa-P 



(Pa,a\p,b) 

a+b { Pa ,a\X \p,b) 
(p a ,a\p,b) 



(65) 



P=Pb 
(66) 



These last two formulae are exact and valid for all values of p a ,p b ,a,b. (They have been 
derived from d6 1 1 under the condition that p a > p b and < b < 1, but as all expressions are 
rational functions of all their arguments, they remain valid everywhere.) 

From ((66) it is possible to show that <P(fi,h) defined by 



(W\exp[{e h D + E)p]\V) 
WW) 



(67) 



satisfies the following equation 

d<P _ b(l+p a {e h -l))+a{l+p b (e h -\)) 



dp 



Pa-Pb 



* + (1 + ^-1))^ 
dph 



This equation can be solved by the method of characteristics, which tells us that the solution is 
of the form 

*b,h) = {l + a ~^ ))b * ((! + P»( e " - ~ !)] 

The fact that <P(0,h) = 1 determines the unknown function & and one gets 



<P(p,h) = 



(p. -/*)(«* -1) 



1 + p a (e h - 1) -exp[/i(e* - \)]{l+p b {e h - 1)) 



exp[^( e "-l)] 



(68) 



(see eq (3.7-3.10) of (9)). 

This expression becomes singular as p — > /Xq with 



Po = 



1 



1 



log 



l+p a (e"-l) 
l+p A («*-l) 



and by analysing the power law singularity one can get the asymptotic expression valid for large 
L 

(W\(e"D + E) L \V) r{a + b + L) {p a - p b ) a + b p Q L -"- b 



{W\V} 



r{a + b) {\ + p a {e h -\)Y {l+ Pb (e l '-l)) b 



(69) 



Second consequence of (SD : 
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Another important consequence which can be obtained by dividing 1 I6U by (62\ is the following 
additivity formula 



{p a ,a\X { X 2 \p b) b) 
{p a ,a\(D + E) L + L '\p b ,b) ' 

(p a -p h ) a+b+L+L ' 



r(L + a + b)r{L' + 1) 
r(L + L' + a + b) 

(fia,a\Ky\p,b) 



Pb<\P\<Pa 



dp 
lilt 



(70) 



(p,\-b\X 2 \p b ,b) 



{pa-p) a+h+L {p-pbY +L ' (p a ,a\(D + E)L\p,b) {p,l-b\(D+E)L'\p b ,b) ' 

which is the same as eq. (65) of 1131 up to the prefactor which was wrong in fT3l and which is 
corrected here. This formula allows one to compute the properties of a lattice of L + L' sites if 
one knows those of two systems of size L and L'. 

Third consequence of (61}: 

Using <66t and i62t one can write an exact recursion for Z/, denned in (2) 



f>M + h J?? ?L A 



L + a + b 



7 , {Pa-Pb)(l+Pb A+l) dZ L 

Zi_l r~, TT ~I7T ( 71 > 

L + a + b apb 



We won't use this recursion relation in this paper, but we believe that it could be an alternative 
starting point to recover the result ((9] and possibly further corrections. 
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B Derivation of (I51I52I53> in the equilibrium case 

Let us consider a site dependent field with small variations 

Zi = hi - h 

around a certain value h. One can then expand Gi defined in { |2l3t in powers of the z,'s 



G L {h ,---h L ) = G L (h, ■■■h) +£z, (m) + - Y^ZiZj {rHnj)c + O (z 3 ) 



where (.) denotes an average in the constant field h. Far from the boundaries i.e. when i S> 1 and 
L — i^> 1, the correlations become translational invariant (because the system is at equilibrium) 



(ni)=g'(h) ; {mnj) c = Cj-i(h) 



(72) 



and 



g "{h) = £ c k (h) 



One can rewrite GlQh,... hi) as 



G L (h u - -h L ) = G L (h, ■ ■■h) + g'(h)'£zi + ^1^ + 1 c k {h) YziZi+k 



i k>\ ;=1 



(73) 

and using the fact (wich follows from J73t by looking at the term proportional to L when all the 
hj's are equal) that 

2 °° 

g(h i )=g(h)+z i g'(h) + ± £ c *(/7) + 0(,- 3 ) 



one gets 

Gi(h u - h L ) - £g(/i,) = G L (h, -h)- Lg{h) 

i 

~L-k 



(74) 



J i=l i=L—k+l 

1 . 



- J> ( («*) - g' (A) J + jE^IW. - c^(fc)J + O (z 



In the large L limit, the correlation functions, near the boundaries, have a limit which is not 
translational invariant 



(nt) - g'(h) -4 c^(h) ; (n L -i) - g'(h) -> 



whereas 



(n,nj) c -cj- t (h) -> *jy(A) ; (n w ji w ) c -cj_,(fc) -> fc*f (A) 
One then should have 



rfA M (A) _ " , eft 



=i 



rfA ri « ht (/i 



V X -1 right/, \ 



(75) 



!=0 
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daf ft (h) 



dh 



^ ght w =I ^« w _ L ct(A) (76) 



J/7 



*'>!+! 



j k>i 

so that using 446b and the fact that c k {h) = C- k (h) 

— ^2— = L M*)-E* c *(*) : J/;2 = L M*) -£*<*(*) 



For large L this becomes 

Gt(Ai, -*/.) " J>(*0 = A left (/ ! )+A n s» 



;">0,/>0 



■E c *(*) 



7 2 +7 2 



1 



-I z i ( <"■■> "S'W) + 5 E**/ [<Wj)c - Cj-i(h)) + O (z 3 ) 
hi 



which can be rewritten, up to terms of third order in the z, 's 

G t {h u -h L )-Y,g(hi) = D ,eft + D^ ht - £ C ,(/ !( ) 

i k>l 

where 



L y k (hi-h i+k ) 2 
i=i z 



D"*=A , « ft (*i)--£ Ct (A 1 ) 



■E^-AO^Ai) 



■J>-Ai)fa-Ai)i$?(Ai) + 0( 



(77) 



and 



■£(/ )i _,-/ !L )af ght (/ !L ) 



All the differences /j, — /i ; which appear in {77} are between nearby sites i,j. Under this form, 
the differences /?, — hj between remote sites do not need to be small. In what follows we will 
assume that {(77) remains true as long as these differences /i, — hj remain small for nearby sites 
(i.e. for \i — j\ <SlX) even if these differences could be large when \i — A). 

Now for a slowly varying field of the form {4}, with L = yX as in {5} one can evaluate the 
different terms using the Euler Mac Laurin formula 



£ g(hi)~X [ y g(H(x))dx 
,ti Jo 



H'(y)g'(H(y))-H'(Q) g '(H(Q)) 
24A 



;A left (//(0))- 



H'lO) 



*>1 



y (h i -h i+k ) 2 
*-> 2 



1 /? 

2A 



[ V ^k 2 c k (H(x))H'(x) 2 dx 
Jo t >, 
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